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ABSTRACT 

We observed the transiting super-Earth exoplanet GJ1214b using Warm Spitzer at 4.5 /xm wavelength 
during a 20-day quasi-continuous sequence in May 2011. The goals of our long observation were to 
accurately define the infrared transit radius of this nearby super-Earth, to search for the secondary 
eclipse, and to search for other transiting planets in the habitable zone of GJ1214. We here report 
results from the transi t monitoring of GJ1214b, including a re-analysis of previous transit observations 
by iPcsert ct al. (1201 It ). In total, we analyse 14 transits of GJ1214b at 4.5 /xm, 3 transits at 3.6 /im, 
and 7 new ground-based transits in the I-|-z band. Our new Spitzer data by the mselves eliminate 
cloudl ess solar composition atmospheres for GJ1214b, and methane-rich models from lHowe fc Burrowsl 
dMl). Using our new Spitzer measurements to anchor the observed transit radii of GJ1214b at long 
wavelengths, and adding new measurements in I-l-z, we evaluate models from'Bennek e fc Seaged (|2012t ) 
and , Howe fc Bur rows (^201^) using a analysis. We find that the best-fit model exhibits an increase in 
transit radius at short wavelengths due to Rayleigh scattering. Pure water atmospheres are also possible. 
However, a flat line (no atmosphere detected) remains among the best of the statistically acceptable 
models, and better than pure water atmospheres. We explore the effect of systematic differ ences among 
results from different observational groups, and we find that the iHowe fc Burrowsl (|2012f) tholin-haze 
model remains the best fit, even when systematic differences among observers are considered. 

Subject headings: Planets and satellites: atmospheres - Techniques: photometric - Infrared: planetary 
systems - Planets and satellites: composition 



1. INTRODUCTION 



The mass and radius of the nearby tran siting super- 
Earth GJ1214b (jCharbonneau et al.l I2009D implv that 



it mu st have a significant atmosphere ()Rogers fc Seageij 
l2010f ). That inference motivated an extensive effort to de- 
tect the atmosphere, by seeking wavelength variations of 
the transit depth. A wide variety of compositions are pos- 
sible f or super-Earth atmospher e s (iMiller-Ricci fc FortnevI 



120101: iBenneke fc Seaged 120121 : iHowe fc Burrowsl l2012fl 



Most 



current observations of the transits ( 


Bean et al.l 20101 


2011 


Crossficld et al.ll2011 


Desert et al. 


l201lHBerta et al.l 


2011 


2012t ide Mooii et al. 


20121) have reiected hydrogen- 



(|201lh concluded in favor of a low molecular weight at- 
mosphere. The infrared (IR) spectral region is partic- 
ularly important for such studies, because strong water 
vapor bands increase the transit depth in the IR signifi- 
cantly as compared to the optical. This is especially true 
for hydrogen-dominated atmospheres, because of the in- 
creased atmospheric scale height. The intrinsically strong 



IR water vapor opacity makes hydrogen-dominated atmo- 
spheres opaque in the IR over several scale heights, in spite 
of their relative paucity of heavy elements. Also, strong 
bands of methane and carbon monoxide fall within the 
Warm Spitzer bandpasses at 3.6- and 4.5 ^m, respectively. 
Considering that GJ1214b's M-dwarf host star is bright 
in the IR, transit observations of GJ1214b using Warm 
Spitzer become particularly relevant to the characteriza- 
tion of its atmosphere. 

In May 2011, we observed the GJ1214 system quasi- 
continuously for 20 days, using Spitzer's IRAC instrument 
at 4.5 /xm. Our investigation had three goals: 1) to im- 
prove the transit parameters of the system and constrain 
the properties of the planet's transmission spectrum, 2) to 
search for the secondary eclipse, and 3) to search for other 
transiting planets in this system, to the outer edge of the 
habitable zone. 

Our observations included minor interruptions for data 
downloads. But, an unplanned 42-hour data loss also oc- 
curred during the 20-day sequence, caused by a combina- 
tion of spacecraft and Deep Space Network (DSN) down- 
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link anomalies. We therefore re-observed GJ1214 for an 
additional 42-hours in November, 2011, using the IRAC 
3.6 /im band. Consequently, we have multiple transits at 
both Warm Spitzer wavelengths, and these data provide a 
particularly powerful constraint on the IR transit depth. 
Moreover, since our 4.5 /.tm transits were observed nearly 
consecutively, we have an excellent basis for evaluating the 
degree to which stellar activity (e.g., star spots) affect the 
inferred transit depth. 

A glob al analysis o f these GJ1214 data is being pub- 
lished bv lGillon et al.l (|2012[ ). including a search for other 
transiting planets in this system. In this paper, we focus 
on the implications of our observations for understanding 
the nature of GJ1214b's atmosphere, and we investigate 
the wavelength-dependent transit radius of GJ1214b in de- 
tail. Our analysis includes the degree to which star spots 
- even those not occulted by the planet - contribute to 
possible bias in the measured radius of the planet. An- 
chored by our improved precision for these infrared tran- 
sits of GJ1214b, we add 7 new ground-based transits in 
the 1+z band, and we re-analyze the totality of pub- 
lished wavelength-dependent transit depths for GJ1214b, 
exploiti ng recent advances in super-Earth model atmo - 
spheres ()Benneke fc Seage"ill2012l : iHowe fc Burrowsl[2012l) . 

Sec. 2 describes the details of our observations, and 
Sec. 3 explains our procedure to extract precise photom- 
etry from the data. In Sec. 4 we fit to the photometry, 
extracting the transit radius for the planet in the warm 
Spitzer and I-|-z bands, and we derive improved system 
parameters. Sec. 5 considers the possible effect of star 
spots on our results. Sec. 6 discusses implications for the 
atmosphere of GJ1214b. 

2. OBSERVATIONS 
2.1. Spitzer 

We observed GJ1214 for 20 consecutive days using 
Warm Spitzer at 4.5 /im, beginning on April 29, 2011 at 03- 
46 UTG. We used subarray mode with an exposure time of 
2 seconds per frame. The observations contain several ~ 3- 
hour interruptions for data download, and one unantici- 
pated 42-hour gap where data were irretrievably lost. The 
data loss occurred because the DSN incurred anomalous 
delays in downloading data at that time. When the Spitzer 
observatory was designed, long observing campaigns of 
exoplanet photometry were not envisioned. Spitzer's on- 
board flight software was designed to automatically delete 
data after a certain period of time, to make room in mem- 
ory for new observations. The DSN anomaly consequently 
caused the onboard software to delete data before down- 
link (the flight software has now been corrected so that 
this will not recur). 

To compensate for the data loss, we were awarded 42 
hours of continuous observations that began on November 
6, 2011, at 11-54 UTG. We elected to acquire these data in 
Spitzer's 3.6 /im bandpass, to complement the 4.5 /im data 
from May 2011. We again used 2-second exposure times 
in subarray mode. In total, our data comprise 791,808 
exposures at 4.5 /im, and 74,624 exposures at 3.6 /im. 

2.2. TRAPPIST 

In order to help define the possible effects of stellar 
activity on the Spitzer transits, to further cross-check 



our analysis versus [Gillon et all (j2012l) . and to add addi- 
tional information relevant to the atmosphere of GJ1214b, 
we observed 7 transits using the TRAPPIST facility 
(iGillon et al.ll2011l:lJehin et al.ll20lH) . over the period 2011 
March 11 - May 18. The T RAPPIST observat ions and 
photometry are described bv iGillon et al.l ()2012[ ). but we 
summarize the data here. The observations were made 
using the 60-cm robotic telescope in a slightly defocused 
mode. An I-l-z filter gave transmission from 750 to 1100 
nm. Differential photometry on the 25-sec exposure im- 
ages was done (by M.G.) using IRAF/DAOPHOT. In our 
analysis, we use t he same version of the photometry as 
IGillon et al.l (|2012[ ). but we perform an independent anal- 
ysis and transit fitting. 

3. SPITZER PHOTOMETRY 

3.1. Aperture Photometry 

Our analysis utilizes the Basic Galibrated Data (BCD) 
files produced by version SIS. 18.0 of the Spitzer pipeline. 
Two dimensional (2D) Gaussian centering prod uces the 
least scatter in our fi nal photometry (Steven son et al.l 
120101: lAgol et al.|[20Tol ). iKnutson et ad (|2Q12[ ) found that 
flux-weighted centering gives superior results for Spitzer 
data in studies of exoplanetary phase curves over long time 
scales. We tried flux-weighted centering for our transit 
analysis, but it did not result in significant improvement 
over our 2D Gaussian centering. 

In the case of the 4.5 /xm photometry we center a circular 
aperture of constant radius on the star. We calculate the 
stellar flux within the aperture, including analytic approx- 
imations for the partial coverage of pixels at the boundary 
of the aperture. We vary the radius of the aperture from 
2.0 to 5.0 pixels, in 0.5-pixel increments, and thereby pro- 
duce seven versions of the photometry at each wavelength. 
After decorrelation (see below), we chose to use an aper- 
ture radius of 2.5 pixels for 4.5 //m, based on the global 
scatter in the decorrelated photometry. 

In the case of the 3.6 /im photometry, we examined con- 
stant radius photometry- for the same seven radii as in 
the 4.5 /im section- and variable radius aperture photom- 
etry. The variable radius photometry improves the preci- 
sion by 41% over constant aperture photometry by varying 
the aperture radius a s a function o f the 'noise-pixel' cal- 
culations per fram e ()Mighelll 120051: iKnutson et all 120121: 
iLewis etall l2013t ). Therefore we adopt the noise pixel 
method for our 3.6 /im photometry, and subsequent intra- 
pixel decorrelation (see below). 

The noise pixel method estimates the effective width of 
the pixel response function, accounting for undersampling, 
by calculating the variance of the flux per frame, weighted 
by the square of the mean flux: 



where pi is the noise pixel parameter for the stellar im- 
age in frame i , Ij is the intensity of pixel j , and the summa- 
tions extend over all pixels wherein the stellar intensity is 

significant. Using the \J~^i as the aperture radius collects 
an optimum amount of light for photometry. The average 
aperture radius from using this formulation on our 3.6 /zm 
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data is 2.60 pixels. 



3.2. Decorrelation 



Upon producing photometry, we immediately see 
Spitzer's well-known intra-pixel sensitivity effect, that 
must be decorrelated and removed from the data. A 
portion of the raw 4.5 fim photometry (before decorrela- 
tion) is illustrated in Figure 1. Because the transits of 
GJ1214b are substantially larger than the intra-pixel sig- 
nature, we first mask off the transits from the first stage of 
decorrelation. Our 4.5 fim data comprise seven distinct re- 
acquisitions of GJ1214, with interruptions for data down- 
load. Therefore our first stage of decorrelation is done 
separately for each of the seven re-acquisitions. That ini- 
tal decorelation fits two polynomials to the data for each 
re-acquisition. (There are seven sections of the data de- 
fined by the data downloads and re-acquisitions.) One 
polynomial fit is applied to those data points wherein the 
image centroid lies at Y-coordinate greater than the center 
of the pixel, and another polynomial fit is applied to those 
data wherein the image lies at Y-coordinate less than the 
center of the pixel. Our rationale for this two-parameter 
fit is based on visual inspection of the photometry, that 
shows different behavior above the center-of-pixel than be- 
low center (see lower panel of Figure 1). The polynomi- 
als are fourth order in Y and second order in X, because 
the photometric variations as a function of Y are more 
pronounced than with X. Also, visual examination of the 
variation with Y indicated that lower order polynomials 
(e.g., quadratic) would not represent the variations opti- 
mally. After finding the best fit polynomials, we divide 
them into the data, including the in-transit data. We only 
used this step to determine initial conditions for an itera- 
tive weighting function method (see below). 

We follow the polynomial decorrelation with a sec- 
ond procedure that is itself a two-pass iterative pro- 
cess. We first remove a preliminary transit model from 
the polynomial-dec orrelated data, and apply a Gaussian - 
weighting function (iBallard et al.|[2010l : iLewis et"al]|2013t 
iKnutson et al.ll2012[ ) to correct the intra-pixel signatures. 

The kernel of the Gaussian weighting uses a cr of 0.005 
pixels in Y, and 0.01 pixels in X for the 4.5 fim data, that 
we determined by trial and error - evaluating the noise 
level of the final decorrelated photometry. A separate 
Gaussian weighting was applied to each section of the data 
between downloads, but the same kernel size was used for 
all data sections at 4.5 /im. 

For the 3.6 fxni data, we applied a noise-pixel, Gaus- 
sian weighting function that uses a variable a in Y, X, 
and /3 (noise-pixel value). This varied the weights in the 
Gaussian kernel to accommodate the necessary number of 
neighboring points that infiuence the strength of the cor- 
relation between c enter posi tion and flux, as discussed in 
ILewis et all ()2013[ ) and .Knutson et al.l (|2012D . 

Because of effects near the pixel boundary, we chose to 
subdivide the 3.6 /im data into 2 sections, encompassing 
the 2 transits in the data. The pixel boundary effects oc- 
cured well out of phase of both transits. After examining 
the global scatter and a vs. N~'^-^ (bin size) slope residu- 
als, we found that decorrelating after subdividing the 3.6 
/xm data set to within ±0.05 of each transit produced the 
best noise levels, and resulted in more conservative un- 



certainty estimates on the transit parameters, by ~ 10%. 
Moreover, this phase range coincides with similar subdivi- 
sions in the 4.5 /im data. 

After dividing by the results of the Gaussian weight- 
ing, we solve for and subtract an improved transit model, 
and again apply a second stage weighting function decor- 
relation. We experimented with a third iteration of this 
process, but it did not produce significant improvement. 
Unlike the first (polynomial) stage of decorrelation, the 
Gaussian weighting function was applied to the totality of 
the data (i.e., all re-acquisitions) using a single Gaussian 
kernel. The initial polynomial stage of the decorrelation 
process may seem unnecessary, since the weighting func- 
tion alone could remove structure in the data on both large 
and small spatial scales. However, we find that the poly- 
nomials speed up the iterative process by providing a fast 
start to the iteration. Our final photometry has a stan- 
dard deviation of 3.7 x 10"'^, which is only 15% greater 
than the photon noise. Moreover, red noise is minimal 
and the precision for binned data improves nearly as the 
square-root of the bin size, as we demonstrate below. 

Figure 2 shows an overview of our 4.5 fim photometry, 
after decorrelation; more detailed depictions are discussed 
and shown below. 

4. MODEL FITTING FOR TRANSIT PARAMETERS 

We analyze all available Sptizer trans its, including a re - 



analysis of the transits reported by D esert et al.l (120111). 
and T RAPPIST transits using the I-l-z filter (jGillon et al I 
I2012D . We use three methodologies to determine the best- 
fit transit parameters for each dataset and wavelength. All 
three methods use only the data within a phase interval 
of ±0.05 around the center of each transit- except for the 
TRA PPIST trans i ts, wh ich used all of the available data 
from IGillon et al.l (|2012l) . 

The first method solves for the best-fit transit parame- 
ters of all transits simultaneously at each wavelength. The 
second method fits each transit individually and indepen- 
dently, then calculates the average of the transit param- 
eters at each wavelength, weighting the individual results 
by the inverse of their variance. The third method phases 
and bins all of the transits at each wavelength into a sin- 
gle transit, and fits to those phased & binned data. All 3 
methods included a total of 3 transits at 3.6 /xm, 14 tran- 
sits at 4.5 /im, and 7 transits in the I+z band. Comparing 
these three methods gives an indication of the consistency 
of our results, and we do indeed find good consistency, as 
noted below. 

We now describe the details of those fitting procedures. 

4.1. Spitzer 4-5 fJ.m Transits 

We use the formulation of lMandel fc Agoll (j2002|) to gen- 
erate transit curves and fit them to the observed data, 
thereby extracting the essential parameters of the tran- 
sit. The 14 Spitzer transits at 4.5 /im comprise our most 
extensive and highest quality data. Our simultaneous 
fit holds the orbital period fixed at the value measured 
bv lBean et aP (l201lD (1.58040481 ± 1.210"^), and uses a 
Levenberg-Marquardt algorithm to minimize the for 
each transit. Although we are minimizing x^, we do not 
accept that particular set of transit parameters as our best- 
fit values. Instead, we explore parameter space using a 
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Markov Chain Monte Carlo (MCMC) method, and also 
using a residual-permutation ('prayer-bead') method. Our 
best-fit values are taken from the medians of the posterior 
distributions generated in this exploration of parameter 
space (see Sec. 4.4). 

Our fit extracts a correction to the transit center time 
(as might be caused by ephemeris error), as well as a/i?*, i, 
Rp/Rif, and a linear limb-darkening coefhcient. We hold 
the quadratic limb darkening c oefficient at zero because 
Bayesian Information Criterion (|Liddlell2007D analysis sup- 
ports a linear law, which gives an adequate account of the 
minimal limb darkening that is characteristic of infrared 
transits. Moreover, our derived linear limb darkening co- 
efhcients at 3.6 and 4.5 microns (cO =0.158 and 0.128 re- 
spectively, see Tabl e 4) are reasonably c onsist ent with the 
values predicted bv IClaret fc BloemenI ()20lH ) (cO= 0.147 
and 0.155) for a model atmosphere having T=3500K (that 
temperature being their closest match to GJ1214). 

Figure 3 zooms in on a portion of the simultaneous fit 
for the first two 4.5 /im transits, and Figure 4 shows the 
standard deviation of residuals (data minus simultaneous 
fit) when binned over different time intervals. Note that 
we find very little red noise in these data, indicating the 
success of our multi-stage intra-pixel decorrelation. Note 
also that our decorrelation process will tend to remove slow 
variations in the stellar brightness. Hence the apparently 
constant flux seen in Figure 2 should not be interpreted 
as evidence for stellar quiesence. (The possible effects of 
stellar activity on our results are discussed in Sec. 5). 

In a ddition, we b i n all 14 4.5 /im transits, including 
that of iDesert et~al] (|201lD . into bins of width 0.001 in 
phase, using a running standard deviation for the weights 
in a weighted mean- phasing them all to a common epoch 
determined from the simultaneous fit. As before, we in- 
clude only data within a phase range ±0.05 of transit 
center, and we fit to the phased & binned transit us- 
ing the same Levenberg-Marquardt algorithm described 
above. Figure 5 shows the resulting fit in comparison to 
the phased & binned data. 

The best-fit Spitzer transit parameters for each 4.5 /im 
transit fitted individually are given in Table 1; the indi- 
vidual 3.6 /im transits (Sec. 4.2) are given in Table 2, and 
the individual TRAPPIST results (Sec. 4.3) are given in 
Table 3. Results using the combining methods are sum- 
marized in Sec. 4.5. 

4.2. Spitzer 3.6 pim Transits 

Our 42-hour 'replacement' observations contain two 
transits at 3.6 /im, one near the beginning of these data 
and one near the end. The photometry for these transits 
was decorrelated using the same methodology described 
above for 4.5 /im. In contrast to the 4.5 /^m case, our 
3.6 /im photometry exhibits noticeable red noise. Fortu- 
nately, this red noise is most significant in the long in- 
terval between the two transits. We limited the effect of 
this red noise by limiting the range of the data included 
in our decorrelations and fits. The omitted data did not 
occur near the transits, nevertheless we tried to develop 
objective criteria for the range of data that were used. 

After comparing the a vs. N~''^-^ and global a, or RMS 
scatter, of various sized data slices, we determined that 
trimming the 3.6 /^m data set into a phase range of ±0.05 



around the center of each transit minimized the RMS scat- 
ter and the residual of the a vs. N~'^-^ slope, which mini- 
mized the red noise for these transits; and was fortuitously 
the same range as the ±0.05 in phase that we adopted for 
our 4.5 /im fits. 

We repea ted all of the method ology described in this sec- 
tion on the iDesert et al.l (j201lh 3.6 /xm transit data, and 
we found that the slope of a vs. N~'^-^ was insensitive to 
the range of data a nalyzed. As a result, we use the entire 
IDesert et all (|2011h 3 .6 /im data set. 

As in the 4.5 /im case, we fit to all 3 of the 3.6 /im transits 
using 3 methods: simultaneously, individually, and phased 
& binned. Figure 6 shows all 3 of the 3.6 fim transits 
phased & binned, overlaid with a best-fit curve. Figure 7 
illustrates the standard deviation of the 3.6 /im residuals 
as a function of bin size for our fits to the simultaneous 
transit parameters. 

4.3. TRAPPIST I+z Transits 

For the TRAPPIST data set, we used aU 7 distinct 
epochs of t h e GJ 1214b transiting system provided by 
iGillon et al.l (|2012t ). Three of these epochs overlapped 
with the Spitzer 4.5 /im data set. iGillon et all (|2012[ ) used 
the I+z filter on the TRAPPIST telescope because it sup- 
plied a near uniform filter profile from 0.7 - 1.0 /im. 

We determined the physical parameters using all three 
methods discussed above: a simultanous fit, individual 
fits, and a fit to phased fc binned data. Simila r to the 
Spitzer data sets, we fit a iMandel fc Aeol (|20Q1 transit 
model using a Levenberg-Marquardt routine. In contrast 
t o the Spitzer data sets, we fit the TRAPPIST data using 
a i Mandel fc Agoll (j2002|) model that included quadratic 
limb darkening, which accounted for the excess stellar 
limb-darkening observed in the shorter wavelength tran- 
sit data. The fit to the phased & binned TRAPPIST data 
is shown as Figure 8. To analyze the quality of the fit. 
Figure 9 shows the a vs. N~^-^ for the simultaneous fit of 
the TRAPPIST data set. 

The results for the planet-to-star radius ratios, and re- 
lated error bars, in the I+z band are plotted on Figure 10 
for comparison to our 13 Spitzer transits at 4.5 /im, and 
listed in Tables 3 & 4. Further information about the data 
re duction process for the TRAPPIST data set is included 
in IGillon et all (I2012D . 

4.4. Errors 

We used two methods to estimate uncertainties for our 
derived transit parameters: MCMC and prayer-bead. In 
both methods the errors - as well as the best-fit values 
- follow from the posterior distributions. We adopted 
the prayer-bead method for our quoted re sults because 
it expl icitly includes the effect of red noise ()Desert et al.l 
l2011bD . Figure 11 shows the distributions for Rp/R^, for 
one of our 4.5 /im individual fits, showing a broader distri- 
bution for the prayer-bead method. 

In implementing the prayer-bead method, we permute 
the residuals only within the adopted phase range of the 
fit (±0.05 in phase). At each permutation, we find the 
best-fit transit parameters using the Levenberg-Marquardt 
method, and we add those best-fit values to the posterior 
distributions of each parameter. We adopt the median of 
the posterior distribution as the best fit value, following 
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iDesert et all ()2011b[ ). 

A potential additional source of error is associated with 
the decorrelation, that is not explicitly propagated into the 
stage of fitting the photometry. However, we are not con- 
cerned about this for two reasons. First, the 4.5 /xm data 
are so extensive, and the intra-pixel effect is so modest at 
that wavelength, that we believe those errors in decorrela- 
tion have negligible effect. Second, our procedure accounts 
for imperfect decorrelation at both Spitzer wavelengths, in 
an implicit fashion. Imperfections in decorrelation create 
red noise, and that red noise contributes to errors on the 
derived transit parameters using the prayer-bead method. 

In addition to random error, systematic differences may 
exist between our results and other investigators. We dis- 
cuss one source of possible systematic difference in Sec. 4.5; 
the implications of these differences for the nature of 
GJ1214b's atmosphere is discussed in Sec. 6. 

4.5. System Parameters 

As noted above, we estimated the system parameters us- 
ing three methods: 1) simultaneous fitting of all transits at 
a given wavelength, 2) averaging system parameters from 
the individual fits to each transit at a given wavelength, 
and 3) fitting to phased & binned combinations of transits 
at a given wavelength. We included the transits observed 
by Desert et al. (2011) in all three methods. We also fit the 
TRAPPIST transits using all three methods (see Table 4). 
Table 1 lists the individual fits to the 4.5 fj,m data; Table 2 
gives the individual fit results at 3.6 /xm and Table 3 lists 
the individual fit results for the TRAPPIST transits in the 
I-l-z band. 

Since our different fitting methods (summarized in Ta- 
ble 4) are simply different ways of accounting for the same 
data, they should give consistent results as far as the best- 
fit parameters are concerned. Beyond best-fit consistency, 
we find that comparison of these methods can provide a 
basis for caution concerning the errors on the derived pa- 
rameters. For example, our 4.5 /im results from the simul- 
taneous fit give Rp/Rf, = 0.11710 ±0.00017, whereas aver- 
aging the individual fits, weighted by the inverse of their 
variances, gives Rp/R^ = 0.11699±0.00026. Affhough the 
best-fit values agree well, the larger error from averaging 
the individual fits may indicate potential variations from 
transit-to-transit . 

Table 5 summarizes our results for the Rp/ R^, parameter 
that potentially reveals information about the atmosphere 
of GJ1214b into one concise Table. For the discussion that 
follows, we adopt the results from the phased & binned 
method, because we feel that the high precision of these 
combined transit curves allows the most reliable solution. 

We explored further comparis on with tw o other pre- 
cise measurement s of Rp/R^: iBean et"al] (|201l[ ) and 
iBerta et al.l ()2012f ). We choose these investigations for 
more in-depth comparison because the former is the high- 
est precision ground-based measurement of GJ1214b, and 
the latter is a precise space-borne measurement. We inves- 
tigated to what extent differences in Rp/R^, arise from the 
different transit solutions, with different values for orbital 
transit parameters such as a/i?* and orbital inclination. 
Because these orbital parameters should not vary with 
wavelength, we force our sol utions to adopt the values as 
derived by , Bean et al.l (|201ll) and .Bcrta et al, (,2012,). Our 



resultant retrievals for Rp/ R^ are included in Table 5. On 
average, constraining o ur orbital parame ters to have the 
values found by either iBean et al.l (|20lH ) or iBerta et all 
(2012) results in decreasing our i?p/i?* value at 4.5 ^m by 
about 0.0010, but with less difference at 3.6 ^m or I-fz. 

Arguably, we should adopt these constrained values as 
our principal result. However, the low limb darkening that 
prevails at Spitzer wavelengths, in combination with the 
high precision we achieve from our large dataset, moti- 
vates us to rely primarily on our own orbital parameters. 
Nevertheless, we explo re the implications of adopting the 
iBean et all (|201lh and lBerta et al.l (|2012D orbital parame- 
ters in Sec. 6. 

5. TRANSIT-TO-TRANSIT VARIABILITY AND STAR SPOTS 

A star spot crossing during transit appears as an 
anomalous spike or bump in the transit light curve (e.g., 
iDeming et al.|[20Tl[) . Our photometry shows no evidence 
that the planet crossed even one significant star spot 
during our thirteen 4.5 /im transits. Nevertheless, spots 
are common on M-dwarf st ars, and uncrossed star spots 
could still affect the transit (jSing et al.ll2009l : iDesert et al.l 
i2011b) . 

TRAPPIST photometery of GJ1214 out of transit, but 
over the same time period as our 4.5 /im Spitzer ob- 
servations, shows essentially no variation in the I-band 
(jGillon et al.ll2012D . to a limit of about 0.2%. Neverthe- 
less, other investigations have found that GJ1214 exhibits 
rotationally-modulated signatures of star spots, so we con- 
sider the potential impact of such variation on our results. 
As we will demonstrate below, even allowin g for more pho- 
tometric variation than lGillon et all (j2012[ ) observed, star 
s pots have negligible effect on our results. 

IBerta et al.r(|201ll ) found that GJ1214 shows a photo- 
metric variation of 1% amplitude (2% peak-to-peak) in the 
I-band (0.715-1.0/im) wi th a rotation period of ~ 53 days , 
based on MEarth data (|Nutzman fc Charbonnea"ij||2008( ). 
spanning three years of observations. Using P HOENIX 
model atmosp heres lAllard fc Hauschild"3 Il995l) , we cal- 
culate that the iBerta et al. ( 20111) amplitude of variation 
could be produced by two star spots, each covering as 
much as 2% of the sky-projected stellar disk, separated in 
longitude by 180°. Based on Doppler imaging studies of 
active dwarf stars (Rice & Strassmeier 1998), we adopt a 
temperature contrast (spot vs. photosphere) AT/T = 0.1, 
thus r^nnt ^ 2700 j-C. Note that this is the same as adopted 



bv IBerta et all (l20TTI) . We use limb darkening coefficients 
for both the photosphere and the star spot, calculated from 
our transit fitting at 4.5 /im (cq = 0.11, ci = 0.0). On this 
basis, we determined that a 1% variation of the stellar 
light curve in the I-band translates to 0.42% variation in 
Spitzer photometry at 4.5 fj,m. 

We developed a numerical tile-thc-star model to calcu- 
late the effect of unocculted star spots on the transit depth. 
We created a synthetic time-series of 2D images of GJ1214, 
at 4.5 /im, and projected two circular star spots on its sur- 
face at the equator, separated by 180° of longitude. This 
arrangement of spots in opposite hemispheres produces an 
appropriate quasi-sinusoidal effect in the total stellar light. 

We accounted for variation in the spots projected area 
as the star rotates, but we ignored the Wilson depres- 
sion effect. The PHOENIX model with T = 3000K, 
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[M/H] = +0.3, log{g) = 5.0 and a = 0.0 represented 
GJ1214; and, the PHOENIX model with T = 27007^, 
[M/H] = +0.3, log{g) = 5.0, a = 0.0 represented the 
star spot. We multiphed Spitzer's 4.5 filter profile by 
the spectral models and integrated over wavelength to cal- 
culate the expected 4.5 ^m flux variations due to stellar 
rotation with the spots fixed in longitude. Both the star 
and spots are affected by li mb da rkening as determined by 
our fitted [Mandel &: Agol ()2002[ ) model parameters. 

It is easy to show that the effect of unocculted star spots 
on the planetary radius derived from the transit is given 
as: 



with e — 



m 



spotless 



spotted 

A, pot 



(1 ~ ^)Iph + ^Ispot 



(2) 



and where ot indexes out of transit, it 
indexes in transit. / indicates the intensity of the stellar 

spot 



disk, Iph being the intensity of the photosphere, and /, 



being the intensity of the spot. 

Using this equation together with flux variations from 
our tile-the-star model, we calculated the potential vari- 
ation in transit depth as a function of time. This quasi- 
sinusoidal variation has an unknown phase because we do 
not know the longitudes of any real star spots on GJ1214. 
Nevertheless, the amplitude of variation in Rp/R^, from 
this model is 9.6 x 10~^, which is negligible compared 
to the observed scatter in our measurements (see Fig- 
ure 11). Moreover, as noted above, the photometric vari- 
ations of GJ1214, observed concurrently with our tran- 
sit data (iGillon et al.l l2012D. were much less than from 
iBerta et al. (2011). We therefore conclude that star spots 
play a negligible role on the observed variations and/or 
possible bias of our inferred radii for GJ1214b. However, 
we cannot exclude the po ssibi lity that incr eased star spots 
during iBean et al.l ( 201lf ) and lBerta et al.l (j201 2) observa- 
tions are responsible for some of the differences between 
our results. 

6. THE ATMOSPHERE OF GJ1214B 

There are numerous transit observations in the 
literature that be ar on the nature of the atmo- 
sphere of GJ1214b ICharbonneau et all 120091 : iBean et al.l 
201 0. 2011: iDesert e t al. 2011; Ber ta et al.l 12011 2012; 
Croll et all 1201 It ISada et al.l l2010l: ICarter et al.l l20Tli 



Kundurthv et al.i201ll:lde Mooii et al.ll2012tlMurgas et all 



2012l : lNarita et alll2012l ). Nevertheless, the improved pre- 
cision we have been able to achieve in the Spitzer bands, 
together with the new TRAPPIST results in the I+z- 
band, and recent advances in modelin g the atmosphere 
of GJ1214b (|Benneke fc Seaged 120121 : iHowe fc Burrows! 
I2012D . motivate us to re- examine the nature of GJ1214b's 
atmosphere. 

In the following sub-sections we review the methodol- 
ogy for comparing observations and models (Sec. 6.1), we 
briefly preview what our Spitzer observations alone can re- 
veal concerning the atmosphere of GJ1214b (Sec. 6.2), and 
we then compare the totality of all published observations 
to existing models, using a analysis (Sec. 6.3). Since 
planetary radii derived by different observational groups 
can differ systematically, we discuss the effect of one par- 
ticular systematic difference in Sec. 6.4, and we summarize 



our conclusions concerning the atmosphere of GJ1214b in 
Sec. 6.5. 

6.1. Comparing Models to Observations 

To compare transmission models for GJ1214b with the 
observations, we need to integrate the models- multiplied 
by the filter profiles- over the observed bandpasses. Let 
Fot(A) be the out of transit flux measured from the star 
as a function of wavelength. Similarly let Fit{X) be the 
in-transit flux as a function of wavelength. Consider the 
simplified case where limb darkening can be neglected (ar- 
guably applicable in the infrared), and include the fact 
that there is a wavelength-dependent observational sensi- 
tivity, S{\). In that case: 

F,t(\) = S{X)h{X){7TRl - 7TRp{Xf), (3) 

where /* (A) is the intensity emergent from the stellar 
atmosphere at wavelength A. The out of transit flux is: 

Fot{X) ^ S{\)F{X)7rRl (4) 

With realistic spectral resolution, the observed quanti- 
ties are integrals over the observational bandpass, and the 
transit depth d is: 



[FotiX) ~ Fu{X)]dX/ / FotiX)dX 



(5) 



Thus: 



d^ j S{X)F{X)Rp{XydX/ J S{X)F{X)RidX. (6) 

When we seek to evaluate a model of the planet's transit 
radius as a function of wavelength (i?p(A)), it is necessary 
to include the wavelength dependence of the stellar inten- 
sity as well as the observational sensitivity. The latter is 
commonly incorporated in numerous studies of both tran- 
sits and secondary eclipses, but the necessity of including 
the stellar intensity is less widely appreciated, especially 
the possible effect of line structure in the stellar spectrum. 
For example, if the stellar spectrum contains water vapor 
absorption that overlaps to some degree with planetary 
water vapor features, then the apparent transit depth will 
be reduced compared to the case where the star is purely 
a continuum source. Stellar intensity weighting is par- 
ticularly important for M-dwarf host stars like GJ1214, 
because their spectrum varies strongly with wavelength in 
the optical and near-IR. Unfortunately, that weighting is 
also considerably uncertain for M-dwarf stars, particularly 
at the very interesting blue wave lengths where the planet 
may exhibit scattering from haze (jBenneke fc Seagerll2012l : 
iHowe fc Burrowsll2012D . 

To compare models of Rp{X) /R^, to observations, we cal- 
culate the value of d using Eq. (6), and infer Rp/R^ as 

^/d. This procedure is valid even at wavelengths where 
appreciable limb darkening prevails. 

Although we use a PHOENIX model atmosphere to per- 
form the weighting over the observed bandpass, this model 
is not ultimately satisfactory for this purpose. For ex- 
ample, in the gree n bandpass (0. 4 6 nm ) where a tran- 
sit was observed bv lde Mooii et al.l (|2012D . the PHOENIX 
model has essentially no flux (many orders of magnitude 
below the peak flux), whereas the real star has suffi- 
cient flux to produce a transit having good signal-to-noise 
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(|de Mooii et al.ll2012[ ). The reason is that the model in- 
cludes only LTE thermal emission from the star, and does 
not incorporate the various emission signatures of mag- 
netic activity. Therefore we can use the PHOENIX model 
only in the red-optical and infrared. Our default procedure 
is to hold the stellar intensity constant for wavelengths 
shortward of 1000 nm (i.e., we set /*(A) = /*(1000) for 
A < 1000 nm), but we verified that using other prescrip- 
tions shortward of this limit (e.g., blackbody spectra) do 
not greatly influence our present results. However, as the 
precision of observations improves, it will eventually be 
necessary to have an accurate spectrum for the host star 
at all wavelengths. 

6.2. Implications from Spitzer 

Prior to an exhaustive analysis of all data versus all 
models, we mention what our new Spitzer data alone 
immediately reveal and/or constrain concerning the at- 
mosphere of GJ1214b. We (BB & SS) generated three 
new model atmospheres for GJ1214b at an equilibrium 
temperature of 546 K, based on the methodology of 
iBenneke fc Seagej (|2012l) . The models are: 1) a H- 
rich solar abundance model, 2) a 'Hot-Halley' composi- 
tion model which begins with solar composition and adds 
minor molecular constitu ents from accreted icy material 
(jBenneke fc Seageij|20l3) . and 3) a pure water vapor at- 
mosphere. Figure 12 shows the result of integrating these 
three models over the Spitzer observational bandpasses, 
and including the stellar intensity using a PHOENIX 
model having Teff/log(g)/[M/H] = 3000 K/5.0/0.3. Based 
on this comparison, the solar composition model is elim- 
inated based on the Spitzer data alone. The water va- 
por model is preferred over the hot Halley model, based 
on the analysis described below. Moreover, we expect 
that methane-rich models having large scale heights (not 
illustrated) will be rejected by the Spitzer data, due to the 
relative lack of an enhanced radius in the 3.6 fiui band - 
that contains the strong 1/3 band of methane. 

6.3. A Analysis 

To gain further qua ntitative insight, we wei ghted 
the three models from IBenneke fc Seageil (|2012[ ) dis- 
cussed above, as well as all of the models given by 
iHowe fc Burro^ (|20T1 over the observed bandpasses of 
all extant transit observations of GJ1214, including the 
TRAPPIST and warm Spitzer data reported here. We fit 
the mo dels to the data us ing a analysis, as described 
below. iBerta et"al] (|2012D applied a similar analysis 
to data from Hubble/WFC3, but not to the complete set 
of data that we do here, and in particular not including 
our new and precise results in the warm Spitzer bands. 
The utility of is well known to be problematic when 
combining data from different observational groups. How- 
ever, x^ is at least an objective way to compare different 
models, and we explicitly consider one cause of observer- 
to-observer systematic differences in Sec. 6.4. 

To evaluate and compare possible models of the plane- 
tary atmosphere, we fit each model to the data by adding 
an adjustable constant (i.e., wavelength-independent off- 
set) to the modeled values of Rp/R^,. We choose the con- 
stant to minimize the of the difference between the ad- 
justed model and the data. Adding this constant is equiv- 



alent to increasing (or decreasing) the size of the opaque 
(i.e., solid) portion of the planet by a small amount. That 
effectively varies the surface gravity of the planet, and 
would strictly speaking be inconsistent with the model 
that is being adjusted. However, the planetary radius at a 
given atmospheric pressure level is not known at a level of 
accuracy comparable to the adjustments we are making. 
Moreover, the requisite adjustments in the model output 
(typically, 0.0005 in Rp/R^), correspond to less than 1% 
differences in surface gravity. We therefore find this proce- 
dure to be a valid and useful tool for testing models versus 
the obs ervations, and we n ote that a similar procedure was 
used bv lBerta et al.l (|2012|). We fit the three models shown 
in Figure 12 ( Benneke fc Seageil 20121), tog ether with all 
of the models from lHowe fc Burrows! ()2012 ). and we cal- 
culate values for each fit. We also fit a flat line to 
the data, i.e., a planetary radius that does not vary with 
wavelength - indicating no signature of the atmosphere. 

Several of the best fitting models are shown in Figure 13, 
compared to the entirety of published radii for this planet. 
In total, there are 97 observations of i?p/i?* versus wave- 
length on Figure 13. The water atmosphere from Fig- 
ure 12, and two models from lHowe fc Burrows! (|2012f ). are 
also included on Figure 13. For the Figure, the models are 
overplotted monochromatically, without integrating over 
the bandpass. However, the x^ values are calculated by 
integrating the model over the bandpass of each obser- 
vation as described above, and adopting the observed er- 
rors from each source. I n the case of multi-band analyses 
(|Kundurthv et al.ll201ll) . we use the total bandpass from 
multiple filters. For the three overplotted models we also 
show the values for the integrals over the Spitzer band- 
passes, as open symbols. 

With 96 degrees of freedom, models can only be rejected 
at the 99.9% confidence level if they have a x^ exceeding 
144.6. Amon g the possible mod e ls. Th e pure water atmo- 
sphere from IBenneke fc Seagerl (|2012f) yields x^ = 142.7 
for 96 degrees of freed om. The hot Halley and s olar com- 



ey a 

position models from I Benneke fc Seagerl (|2012t ) have x' 
values of 167.9 and 1054.7 respectively, and are unlikely 
descriptions of GJ1214b's atmosphere. The solar compo- 
sition model in particular is strongly rejected unless high 
clouds are included (see below). In this respect, we note 
that some discussion of a solar composition atmosphere 
has occurred with respect to transit o bservations near 2 /xm 
(jCroll et al.ll2011l: iBean et al.ll2011l) . We here emphasize 
that a cloudless solar composition model is incompatible 
with the Spitzer data alone, as well as with the totality 
of the observations over all wavelengths. The issue of the 
transit depth near 2 /xm - while an important datum - is 
not crucial to rejecting a solar composition atmosphere. 

Our improvement in the observed Spitzer precision at 
3.6 /im - overlapping the strong 1^3 band of methane, 
pr ompts us to invest i gate th e methane-composition models 
of iHowe fc Burrowi (|2012[ ). All of their methane compo- 
sitions have x^ above 184.6 (the value for 1% methane). 
However, a solar composition model having high clouds 
(down to pressures of 0.1 mbars) produces an accept- 
able x^ — 145.5. This model is plotted as the red line 
on Figure 13; it predicts an increase in radius in the 
difficult-to-observe regions near 2.7 and 3.3 /im. The best 
fit to the data is the yellow line on Figure 13, having 
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= 119.8. This model, from iHowe fc Burrowsl (|2012[ ). 
contains a dense haze of smah (0.1 fim) thohn particles, 
extendin g to very high altitudes (1 /ibar). It was disfa- 
vored bv lHowe fc Burrow s (2012), but our methodology is 
different in that we incorporate a marginally additive con- 
stant when fitting the models. However, we point out that 
also among the acceptable models is a flat line (not illus- 
trated on Figure 13). This null hypothesis (no atmosphere 
detected) yields = 137.0. 

6.4. Systematic Differences Between Observers 

One further check on the acceptable models is to fo- 
cus on the difference between our precise Spitzer radii and 
very precis e radii deri ved at other wavel engths, i.e., by 
iBean et al.l ()2011f ) and iBerta et all (|2012[ ). As noted at 
the end of Sec. 4.5, we fix the variable orb ital param- 
eters (a /i? ^ and i) at t he va lues derived by iBean et al.l 
(|201lh and lBerta et al.l ()2012[ ) and derive a lower value of 
Rp/R* at 4.5 ^m by about 0.001 (Table 5). Using that 
alternative value in our x^ analysis is one way of evaluat- 
ing the possible effect of observer-to-observer differences 
in radii. This procedure increases the x^ values of most 
models, but by varying amounts. Interestingly, the two 
least affected models are the flat line, whose x^ increases 
by = -1-1.6, and the thohn-haz e model. The water 
model from Bc nneke fc Seagerl ( 20121 ) has Sx^ = -1-8.4, and 
their hot Halley model has S\'^ = -1-19.8. The s cattering 
tholin-haze model from iHowe fc Burrowi ()2012|) remains 
as the lowest absolute x^, and has Sx^ = —2.3, i.e. it be- 
comes more likely, not less likely. We conclude that, even 
when systematic differences among observational groups 
are considered, a scattering atmosphere is currently the 
best estimate for GJ1214b. However, we are unable to 
explore the vast phase space of other possible sources of 



systematic error. Hence we also conclude that the null 
hypothesis (no atmosphere detected) remains among the 
most favored models, especially when systematic errors are 
considered. 



6.5. Summary of Implications for the Atmosphere of 

GJ1214b 

We have obtained new radii for GJ1214b in the I-l-z 
band using TRAPPIST, and very precise radii at 3.6- 
and 4.5 /im using Warm Spitzer in a long series of new 
observations. Our x^ analysis indicates that the best-fit 
model for the atmosphere of GJ1214b contains a haze of 
small particles extending to high altitudes, although pure 
water vapor models remain a possibility. However, a flat 
line is among the best-fitting models, particularly when 
observer-to-observer systematic differe nces are cons i dered . 
Therefore we extend the conclusion of iBerta et al.l (j2012|) 
concerning the flatness of the transmisison spectrum from 
1.1- to 1.7 /im to include our new high-precision Spitzer 
measurements at 3.6- and 4.5 /zm. The atmosphere of 
GJ1214b is not unequivocally detected at this point in 
time. 
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Fig. 1. — The top panel illustrates raw data for the first 2.9-day portion of our 4.5 fim photometry, containing two transits. Each data point 
is aperture photometry before decorrelation, binned by every 10 data points. The middle panel shows the Y-center position as a function 
time; that which we correlate with the flux to remove the intra-pixel effect, and related systematic noise sources. The lower panel shows the 
same photometric values, except with the transit regions excluded, plotted versus the Y-pixel position of the stellar image. Note the different 
spatial dependence of the photometry on each side of pixel center. The red curves are the 4th order polynomial fits that we use to initiate 
the decorrelation process (see text). 
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Fig. 2. — Overview of our 4.5 photometry, after decorrelation and binned in 100 two-second exposures per plotted point . The dashed 
hne shows the transit depth that corresponds to one Earth radius. The 13 transits of GJ1214b are apparent. See lGillon et al] I I2012I) for an 
analysis of other possible transiting planets in this system. 
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Fig. 3. — Expansion of a two-transit portion of the 4.5 fim photometry (Figure 2), with the best-fit transit curves overlaid, binned by every 
10 data points. The red lines in the lower panel show a ilc envelope. 



12 



Fraine et al. 




Fig. 4. — Standard deviation of the residuals (data minus the 'simultaneous' fit, in units of the stellar flux) for all 14 transits at 4.5 /im, 
versus bin size. The red lines in the lower panel show a itlcr envelope. 
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Fig. 5. — Phased Sz binned transit of GJ1214b at 4.5 /xm from our 13 Spitzer transits, plus the transit observed by I Desert et al.l l|2011l ). 
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Fig. 6. — Phased Sz binned transit of GJ1214b at 3.6 /xm from our two Spitzer transits, plus the transit observed bv lDesert et al.l l|201lD . 
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Fig. 7 . — Standard devia tion of the residuals (data minus the 'simultaneous' fit, in units of the stellar flux) for our two transits at 3.6 ^m 
plus the lDesert et al.l iWHt) transit, versus bin size. The extra scatter in this figure, compared to Figure 4 above, is related to the number of 
data points binned overall. There are almost an order of magnitude less points at 3.6 fim than 4.5 /im. 
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Fig. 8.— Phased & binned transit of GJ1214b in the I+z band, from TRAPPIST, observed bv lGillon et al.l l|2012l l. 
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Fig. 9. — Standard deviation of the residuals (data minus the 'simultaneous' fit, in units of the stellar flux) for our seven transits in the 
I+z band on TRAPPIST. The extra scatter in this figure, compared to Figure 4 & 7 above, is related to the number of data points binned 
overall. We were only able to bin out to 300 data points ber bin because the most extensive TRAPPIST transit comprised 300 data points. 
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Fig. 10. — Ratio of planet to stellar radius versus epoch for our 4.5 /^m transits, shown in comparison to transits from TRAPPIST. The 
red dotted line shows the median value of our 4.5 fira transits. 
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Fig. 11. — Histograms (posterior distributions) of Rp/R^ for the individual fit to the first of our 4.5 transits, showing results from both 
the MCMC and prayer-bead methods; with dashed lines for the Icr threshold of each distribution. 
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Fig. 12. — Comparison between our Spitzer results and three recent models for the atmosphere of GJ1214b, from lBenneke &: Seaged II2012I) : 
a cloudless solar model, a comet-like model, and a water vapor model. 
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Fig. 13. — Comparison between all published observations of GJ1214b, including our Spitzer and TRAPPIST transits, and recent models. 
Three models having the lowest are illustrated: a water vapor model from Bennckc & Scager (2012) (blue curve), a solar abundance model 
with very high cloud tops (at 0.1 mbars, red curve), and a model with 1-percent wa ter vapor, and a thick t holin haze of 0.1/im particles 
extending to very high altitude (1 /xbar, yellow curve). The latter two models are from I Howe fc Burrowj ||2012I 1. 
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Table 1 

Results of fitting to 13 individual transits of GJ1214b at 4.5 /im, plus the transit observed by IDesert et al] 

(120111 '). 



Epoch 


-: 


2450000(5 J_Dtds) 




i 








a/R, 






Rp/R* 






Co 






211 


5314, 


.214114 


± 0.000085 


88.81 


± 


0, 


.39 


14. 


,92 ±0, 


.40 


0. 


,11726 


±0. 


,00105 


0. 


133 


± 


0, 


,074 


443 


5680, 


.867846 


± 0.000059 


89.30 


± 


0, 


.33 


15. 


,51 ±0, 


.28 


0. 


,11735 


±0. 


,00064 


0. 


,097 


± 


0, 


,046 


444 


5682, 


.448272 


± 0.000067 


88.89 


± 


0, 


.37 


14. 


,98 ±0, 


.39 


0. 


,11692 


±0. 


,00077 


0. 


,208 


± 


0, 


,047 


445 


5684, 


.028698 


± 0.000056 


89.88 


± 


0, 


.36 


15. 


,72 ±0, 


.27 


0. 


,11681 


±0. 


,00057 


0. 


,126 


± 


0, 


.041 


446 


5685, 


.609139 


± 0.000048 


89.11 


± 


0, 


.23 


15. 


,24 ± 0, 


.24 


0. 


,11625 


±0. 


,00034 


0. 


,118 


± 


0, 


.029 


447 


5687, 


.189632 


± 0.000065 


88.23 


± 


0, 


.26 


14. 


,11 ±0, 


.41 


0. 


,11796 


±0. 


,00054 


0. 


,059 


± 


0, 


.051 


448 


5688, 


.769880 


± 0.000094 


89.42 


± 


0, 


.36 


15. 


.38 ± 0, 


.30 


0. 


,11586 


±0. 


,00100 


0. 


,170 


± 


0, 


.067 


449 


5690, 


.350385 


± 0.000052 


88.34 


± 


0, 


.23 


14. 


.33 ± 0, 


.36 


0. 


,11803 


±0. 


,00052 


0. 


,098 


± 


0, 


.050 


450 


5691, 


.930794 


± 0.000067 


88.22 


± 


0, 


.25 


14. 


.24 ± 0, 


.43 


0. 


,11865 


±0. 


,00049 


0. 


,033 


± 


0, 


.033 


451 


5693, 


.511181 


± 0.000058 


88.81 


± 


0, 


.25 


14. 


,95 ± 0, 


.29 


0. 


,11737 


±0. 


,00053 


0. 


,130 


± 


0, 


.046 


453 


5696, 


.671937 


± 0.000073 


88.68 


± 


0, 


.37 


14. 


,75 ± 0, 


.41 


0. 


,11606 


±0. 


,00075 


0. 


,189 


± 


0, 


.047 


454 


5698, 


.252258 


± 0.000047 


88.93 


± 


0, 


.28 


15. 


,14 ±0, 


.30 


0. 


,11567 


±0. 


,00047 


0. 


,056 


± 


0, 


.029 


455 


5699, 


.832550 


± 0.000068 


88.24 


± 


0, 


.28 


14. 


.17 ±0, 


.43 


0. 


,11802 


±0. 


,00062 


0. 


,070 


± 


0, 


.058 


456 


5701, 


.413161 


± 0.000059 


89.61 


± 


0, 


.30 


15. 


.41 ± 0, 


.23 


0. 


,11549 


±0. 


,00066 


0. 


,205 


± 


0, 


.042 



Table 2 

Results of fitting to two individual transits of GJ 1214b at 3.6 ^m, plus the transit observed by IDesert et alI 

(Enni). 



Epoch 


Tc - 2450000 (BJDtds) 


i 


a/ Rt 


Rp/ R„ 


Co 


210 


5312.633724 ± 0.000072 


89.998 ± 0.324 


15.659 ±0.160 


0.11665 ±0.00082 


0.200 ±0.061 


564 


5872.096721 ± 0.000094 


89.144 ± 0.383 


15.331 ±0.354 


0.11619 ±0.00121 


0.114 ±0.065 


565 


5873.677187 ± 0.000062 


89.221 ±0.296 


15.418 ±0.281 


0.11584 ±0.00066 


0.135 ±0.049 



Table 3 

Results of fitting to 7 individual transits of GJ1214b in the I+z-band, observed by TRAPPIST (jGiLLON et AL.|[2bl2l 'l 



Epoch Tc ~245000Q{BJDtdb) i a/R* Rp/R, cp a 

412 5631.875656 ±0.000092 88.035 ± 0.365 13.743 ± 0.443 0.11680 ± 0.00135 0.237 ± 0.094 0.073 ±0.177 

424 5650.840034 ± 0.000152 87.805 ± 0.384 13.215 ± 0.515 0.12389 ± 0.00216 0.331 ±0.176 0.110 ±0.220 

436 5669.805222 ± 0.000110 88.990 ± 0.436 14.945 ± 0.562 0.11389 ± 0.00275 0.422 ±0.170 0.267 ± 0.345 

441 5677.706878 ±0.000170 88.965 ± 0.412 14.923 ± 0.550 0.11573 ± 0.00331 0.428 ±0.198 0.008 ± 0.370 

443 5680.867699 ± 0.000131 88.336 ± 0.397 14.049 ± 0.486 0.12138 ± 0.00182 0.012 ±0.153 0.602 ± 0.298 

453 5696.672075 ± 0.000135 89.873 ± 0.671 14.932 ±0.741 0.10802 ± 0.00296 0.646 ±0.142 0.141 ± 0.242 

455 5699.832535 ± 0.000143 88.360 ± 0.460 14.516 ± 0.641 0.11877 ± 0.00228 0.343 ±0.138 0.005 ± 0.268 



Table 4 

Results of pitting to composite transits, i.e.: simultaneous pits, averaging individual transits, and fitting to phased & binned transits. 



Method / Channel 


'Ic - 2450000(B J Dtdb) 


I 


a/R, 


Hp/ li» 




Cl 



Sunultaneous 3.6 ^jm 4966.524895 ± 0.000047 89.566 ±0.139 15.607 ±0.104 0.11607 ± 0.00030 0.156 ± 0.022 fixed at 0.0 

Sunultancous 4.5 4966.524918 ± 0.000030 88.794 ±0.061 15.049 ± 0.095 0.11710 ± 0.00017 0.110 ±0.014 fixed at 0.0 

Simultaneous TRAPPIST I+z 4966.524881 ± 0.000063 88.240 ± 0.230 14.011 ±0.288 0.11873 ± 0.00111 0.320 ± 0.097 0.187 ±0.176 

Averaged mdividual 3.6 1966.521898 ± 0.000076 89.469 ± 0.221 15.564 ± 0.080 0.11616 ± 0.00019 0.148 ±0.021 fixed at 0.0 

Averaged mdividual 4.5 4966.524925 ± 0.000027 88.803 ±0.127 15.090 ±0.137 0.11699 ± 0.00026 0.110 ±0.015 fixed at 0.0 

Averaged individual TRAPPIST I±z 4966.524968 ± 0.000077 88.437 ±0.057 14.207 ± 0.239 0.11793 ± 0.00182 0.324 ± 0.068 0.146 ± 0.073 

Phased & binned 3.6^111 4966.524156 ± 0.000043 89.469 ± 0.143 15.547 ± 0.091 0.11602 ± 0.00055 0.158 ± 0.035 fixed at 0.0 

Phased & binned 4.5 Atni 4966.524297 ± 0.000016 88.439 ± 0.061 14.442 ± 0.093 0.11709 ± 0.00022 0.128 ±0.018 fixed at 0.0 

Phased & binned TRAPPIST I±z 4966.524122 ± 0.000079 88.314 ±0.177 14.053 ± 0.259 0.11803 ± 0.00079 0.263 ±0.107 0.316 ±0.149 
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Table 5 

Summary of results for Rp/R, derived at different wavelengths using different methods. The phased & binned 

TRANSIT RESULTS FOR Rp / Rt ARE GIVEN FOR SEVERAL CASES: OUR BEST-FIT ORBITAL P ARAMETERS ja/R*, i ) , AS WELL AS 
Rp/R, WHEN WE FORCE THE ORBITAL PARAMETERS TO HAV E THE VALUES DETERMINED By IBeAN ET AL.I (|201lh (14.97, 88.94) 

AND IBerta et AL.I (j2012|) (15.31, 89.30), respectively. 



Method Rp/R, 



Spitzer 4.5 


Simultaneous 


0.1171 


± 


0.0002 


Averaged Individual 


0.1170 


± 


0.0003 


Phased & binned Best-fit orbital 


0.1171 


± 


0.0002 


Bean orbital 


0.1161 


± 


0.0003 


Berta orbital 


0.1160 


± 


0.0003 


Spitzer 3.6 


Simultaneous 


0.1161 


± 


0.0003 


Averaged Individual 


0.1162 


± 


0.0001 


Phased & binned 








Best-fit orbital 


0.1160 


± 


0.0006 


Bean orbital 


0.1163 


± 


0.0002 


Berta orbital 


0.1161 


± 


0.0002 


TRAPPIST l+z band 


Simultaneous 


0.1187 


± 


0.0011 


Averaged Individual 


0.1179 


± 


0.0018 


Phased & binned 








Best-fit orbital 


0.1180 


± 


0.0005 


Bean orbital 


0.1177 


± 


0.0005 


Berta orbital 


0.1172 


± 


0.0005 



